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We study dynamical mass generation in pure Yang-Mills theory and report on a 
recently developed ansatz that exactly solves the tower of Dyson- Schwinger equa- 
tions in Landau gauge at low Euclidean momentum, featuring enhanced gluon-gluon 
vertices, a finite ghost-gluon vertex in agreement with an old argument of Taylor, 
and an IR suppressed gluon propagator. This ansatz reinforces arguments in favor 
of the concept of a gluon mass gap at low momentum (although the minimum of 
the gluon's dispersion relation is not at zero momentum). As an application, we 
have computed the spectrum of oddballs, three-gluon glueballs with negative parity 
and C-parity. The three body problem is variationally solved employing the color 
density-density interaction of Coulomb gauge QCD with a static Cornell potential. 
Like their even glueball counterparts, oddballs fall on Regge trajectories with sim- 
ilar slope to the pomeron. However their intercept at t = is smaller than the 
lo Regge trajectory and therefore the odderon may only be visible in experimental 
searches (for example at BNL) with higher —t than conducted to date at DESY. 
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In a free Yang-Mills theory gauge boson mass terms break gauge invariance and cannot 
be part of the Lagrangian. Moreover, if self-interactions are incorporated perturbatively, 
to any finite order in perturbation theory the radiative corrections do not generate mass 
terms. However non-perturbativc effects can potentially provide dynamical suppression 
of the gauge boson propagator via a dressing function 

(D 

If this function Z is proportional to p 2 in the limit p — > then the complete propagator 
behaves as that of a Yukawa exchange of finite mass. This behavior was conjectured by 
Cornwall pQ. Dynamical suppression in gluodynamics, unlike the electroweak Lagrangian, 
is not facilitated by coupling to an external Anderson-Higgs field, but can only be an 
intrinsic phenomenon due to the non-linearity and strength of the self-coupling. 

The fact that the positions of resonances in the hadron spectrum seem to follow a 
discrete pattern is a clear indication that the gluon degrees of freedom entail a mass gap 
of order 800 MeV. 

This has been previously incorporated in phenomenological models using a constituent 
gluon mass E] , a waveguide cut-off frequency such as the bag model 0] , or a quantum 
string excitation in the flux tube model [5]. 

In the last years the suppression of the gluon propagator has been put on firm theo- 
retical footing in the framework of the Yang-Mills Dyson-Schwinger equations [5[ [7J |S1 El , 
and has also received support from numerical lattice calculations ^01 ID] Q2| • 

The studies in El E3 , employing sophisticated ansatze for the interaction vertices 
were able to provide analytical and numerical solutions for the infrared behavior of the 
Yang-Mills ghost and gluon propagators. The novelty was that the dressing function 
for the gluon propagator in Landau gauge was more suppressed than p 2 7 leading to an 
infinite gluon mass at zero momentum. Conversely, the dressing function for the ghost 
propagator turned out to be infrared singular instead of suppressed. In terms of a real 
constant k E (0.5,0.6) (the calculation of [Hj yields 0.596), they read 

lim Z{p 2 ) = ct-(p 2 ) 2K lim G{p 2 ) = ct-(p 2 )- K . (2) 

In UJ further progress has been made by observing that one can make a self-consistent 
ansatz that simultaneously solves all of the Yang-Mills wave equations in the far infrared. 
The result for the dressing of the Green's function with m gluon legs and 2n ghost legs, 
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A[m, 2n], is 

lim A[m, 2n](p 2 ) = ct ■ (p 2 )(™-™)« . (3) 

allp^O 

By "dressing" we understand that the canonical dimension of the Green's function is 
still multiplying the result, and by "all" that all the momenta entering the scattering 
amplitude need to tend to zero simultaneously. A proof employing the skeleton expansion 
has been given in [Jj and [HJ- In addition to the aforementioned gluon suppression and 
ghost enhancement, one finds from eq. © that the ghost-gluon vertex (n = 1, m = 1) is 
IR finite in agreement with an argument of Taylor |13| and that the three and four-gluon 
vertices are enhanced by (p 2 )~ 3fi and (p 2 )~ 4K dressings respectively. These relations 
between the dressings force cancellations such that all the effective charges defined from 
the primitively divergent vertices in the Lagrangian are finite in the IR, and a s presents 
an IR fixed point, which theoretically supports approaches based on conformal theory 
[H] and is compatible with analytical perturbation theory [151 Ej • 

The demonstration that eq. © indeed solves the tower of DSE's has been given in 
terms of the skeleton expansion. Here we give an argument that is independent of this 
expansion, sidestepping convergence issues. 

We will examine in turn all the contributions to the RHS of the Dyson Schwinger 
equation for A[m, 2n], assume eq. (0 for the coupled-channel Green's functions appearing 
there, and show that the leading IR divergence yields again our counting rule eq. J2J). 

There are several formulations of the Dyson-Schwinger equations, we employ a very 
intuitive one based on a complete resummation of the perturbative series (without the in- 
termediate skeleton expansion), that can conceivably be generalized by analytical contin- 
uation, and one should remember that the DSE's can be directly formulated as identities 
in the path-integral formalism. 

We will draw a box for the Green's function with m gluon and n ghost legs labeled 
by the indices m, 2n. Its DSE reads 

(4) 



m,2n 



Tree + (2a) + (26) + (3)<S(m) 



The tree-level diagrams are the lowest order in perturbation theory and have their 
canonical dimension. All others involve loop integrals and can feature an anomalous 
infrared exponent. Diagrams of type (2a) have at least one external gluon leg coupled to 
a vertex that is not attached to any other external leg. These diagrams can be resummcd 
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to all orders to one of the three following classes (we write below the corresponding 
infrared exponents). All propagators are understood as fully dressed and carrying their 
dimensions given by eq. J3J). 





m — 1, 

2(71+1) 




to + 1, 2n 




m + 2, 2n 



n(n — to + 3) re(n — to + 4) 



The ghost-triangle type diagram on the left is the IR-dominant and yields the behavior 
of the A[m, 2n] function in eq. (0J. Diagrams marked (2b) (see sketch below, left diagram) 
arise due to the four-gluon coupling, that allows two external gluon legs linked together. 
Finally we have diagrams of type (3), that we define as those with no external gluon 
legs (to — 0) that require separate treatment. These are the only loop diagrams in the 
Dyson-Schwinger equations for multighost scattering kernels. Since there is only one 
bare ghost-gluon vertex in the Lagrangian these can all be resummed to the diagram on 
the right (that is leading and again yields our IR counting): 





to + 1, 2n 



n(n - 



TO 



4) n(n — to) 



This exhausts the possible types of proper diagrams that can be resummed from the 
perturbative Yang-Mills series. 

Once we have established that the ansatz eq. Q is a valid solution of the Dyson- 
Schwinger equations in Landau gauge, we turn to phenomenological applications. These 
are based in the sister Coulomb gauge, also transverse but with a Hamiltonian formulated 
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at equal time. The obvious phenomenological consequence of the gluon mass gap is that 
hybrid mesons (qqg), glueballs, and other states whose leading wavefunction in Fock 
space involve constituent glue, will be rather massive, in agreement with estimates in 
lattice gauge theory ^| ^] . We complement lattice studies in that our semianalytical 
methods in the continuum can readily be used to study higher spin states, not easily 
accessible in the lattice. Since glueballs do not have exotic quantum numbers they are 
the most difficult unconventional states to be disentangled from the meson spectrum. 
However they leave some trace in Regge phenomenology. There is no doubt among 
theorists that Regge behavior is present in QCD, but the BFKL, BKP equations are 
valid only for s » —t >> Aqcd, whereas the experimental evidence so far has been 
gathered at small or zero —t (mostly in the form of total cross sections), therefore non- 
perturbative approaches are necessary to study Regge trajectories at low t. The time- 
honored approach to obtain the intercept of a Regge trajectory at t = M 2 = has been 
to study the resonances on the positive (M 2 , J) quadrant of a Chew-Frautschi plot. 

We address these states by means of a model Hamiltonian employing the Coulomb 
gauge formulation and fields, where the kernel is approximated by a charge-density to 
charge-density interaction with a Cornell potential. In a first approximation one neglects 
(suppressed) transverse gluon exchange, that can later be incorporated in perturbation 
theory (a model computation along these lines is available JH1)- In the gluon sector this 
effective QCD Hamiltonian is 

H 9 eff = Tr J dx [IF (x) • IF (x) + (x) ■ (x)] 

- \ J rfxdy^(x)F(x,y)p«(y) , (5) 

with color charge density p^(x) = / o6c A 6 (x)-II c (x), gauge fields A a , conjugate momenta 
II a = — E a and Abelian components BJ4 = V x A a , for a = 1,2, ...8. The Fourier 
transform to momentum space yields 

A°(x) = /(^^[a Q (q) + a»t ( - q)]e -, (6) 
n a (x) = ^|^^[a a (q)-a a t(- q) ] e ^, (7) 

satisfying the Coulomb gauge transverse condition, 
q-a a (q) = 

(— l) AI <7 A1 a a : A1 (q) = 0. Here a°(q) (fi = 0,±1) are the bare gluon Fock operators from 
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which, by a Bogoliubov-Valatin canonical transformation, the dressed gluon or quasi- 
particle operators, a°(q) = cosh 0(g) a„(q) + sinh0(g) a^(— q), emerge. This similarity 
transformation is a hyperbolic rotation similar to the BCS fermion treatment. These op- 
erators excite constituent gluon quasiparticlcs from the BCS vacuum, |f2 >bcs> and sat- 
isfy the transverse commutation relations, [a"(q), a^(q')] = <5 a h(27r) 3 <5 3 (q — q')£> /J!V (q), 
with D tiL ,(q) = (s^u — (— l^ ^pF^ Y Finally, the quasiparticle or gluon self-energy, 
uj(q) = qe~ 2 ®^ q \ satisfies a gap equation j20 : . A more sophisticated treatment is also 
available in the literature |21j . 

In this report we concentrate on glueballs (hybrid mesons can be found in |2S1)- Two- 
gluon glueballs |2U have been variationally computed in the Russell-Saunders (or L — S) 
scheme. The ground state is the L = 0, S = scalar J PC = ++ that we employ 
to fix the cutoff regulating the divergence in the gluon self-energy, to agree with the 
calculation of Morningstar and Peardon J7j. The rest of the parameters being fixed 
from the quark sector, all other glueball masses are predictions of the approach. If 
we plot the spectrum in a Chew-Frautschi plot to reveal Regge trajectories, this state 
cannot be on the leading (left-most) trajectory since it has total J = 0. Its orbital 
excitations with S = 0, L = J = 2,4... fall on a linear Regge trajectory with slope close 
to the pomeron. The second state is the 2 ++ in which both spins are aligned, and this 
one does fall on the pomeron Regge trajectory as observed first by Simonov (0 and 
references therein). This is natural since the intercept of a Regge trajectory with a line 
of integer J usually (but not always) produces a resonance in the spectrum. However, 
the pomeron having a small slope (0.2-0.3 by modern fits [221 ) cannot be represented by 
ordinary mesons. Casimir color scaling between the gluon-gluon and the quark-antiquark 
interaction provide for an enhanced string tension in the glueball system that supports 
the pomeron-glueball conjecture (the slope a'(t) oc 1/er). Turning now to the three body 
problem, the three-gluon variational wavefunction is 



\^ JPC ) = J dqidq 2 dq 3 (5(qi + q 2 + q 3 ) (8) 

^™ (qi, is. q 3 )c Qbc < (qiX teX MI^bcs , 



with summation over repeated indices. The color tensor C abc is either totally antisym- 
metric f abc (for C = 1) or symmetric d abc (for C = -1). Boson statistics thus requires the 
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C = -1 oddballs to have a symmetric space-spin wavefunction taken here to have form 

^w^toi'^qs) = [ci2/(gi,g2) + 023/(02,33) + 

ci 3 /(gi,93)][^ A (qi) + iL A (q 2 ) + y£ (q 3 )] , (9) 

C12 = (l/iil/x 2 |s/^)(s/i s l//3|S/i)(£AS/x|J.M) i (10) 

sufficient to analyze the lightest states. The other two coefficients in Eq. @ can be 
obtained by permuting the indices in Eq. IjlOjl . Several forms for the variational radial 
wavefunction, f(q,q'), involving two variational parameters, /3 and /?', were investigated 
including a separable form, f(q)f(q'), which only involves one parameter. From previous 
experience \'2'S\ . reliable, accurate variational solutions can be obtained if these functions 
have a bell-shaped form with scalable variational parameters. 

The -I — oddballs are not reported in this paper since they likely require the use 
of three body forces |25|. The ground state of the three-gluon system is calculated in 
our approach to be a h glueball (because of color, an annihilation diagram in 
cannot contribute in the odd-C sector), but this is heavier than the corresponding p- 
wave excitation in the two-body system. It is interesting to note also that these 
glueballs appear in the spectrum as apparent r\ c excitations which have been reexamined 
in recent B- meson factory experiments, leading to a revised rj c (2S) mass. Finally we 

present our (PC = ) spectrum in the figure B For comparison we also show the good 

agreement with lattice computations. In addition we give an existing constituent gluon 
model calculation [2] and another by us. The constituent model treatment of does 
not report on higher spin states and we leave it out. One can immediately observe that 
the leading oddball Regge trajectory, in analogy with the pomeron, the odderon, passes 
by the 3 trajectory (three gluons, all spins aligned) and not by the 1~~ as erroneously 
assumed in the literature. We look forward to further lattice computations for higher 
spin states to confirm this point. 

Also shown in the figure is the conventional lu Regge trajectory. As can be seen 
and was to be expected due to the smaller color factor, it has much larger slope but 
higher intercept. At t — it will clearly dominate scattering amplitudes and the odderon 
should not be visible, as is the case in low —t searches at DESY - Although in 
our best variational model the odderon will not be observable (by the time the two 
trajectories intersect, t — —2.3GeV 2 and J < —1, a linear extrapolation from timelike 
t is no more valid), in our approach, we can not categorically preclude future searches, 
either at DESY or BNL, discovering an odderon signature. Such searches, however, 
should also be conducted at suffuciently higher -t (greater than .5 GeV 2 ). Nevertheless 
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we remain pessimistic about the prospects for observing the odderon especially since a 
recent extension by Kaidalov and Simonov of their approach to higher spin states concurs 
with our findings |2*9] . 

In concluding we note this is a period of timely opportunities for exciting phenomenol- 
ogy related to dynamical gluon mass generation and propagator suppression. For exam- 
ple, the recent CLEO measurement [2D] of the l" 1 meson mass, almost degenerate with 
the 1 ++ , shows that the hyperfine interaction is quite short range. This has traditionally 
been a conundrum of the quark model. If a long range Coulomb form is supposed for the 
physical gluon exchange and AV(r) cx —S(r) the hyperfine potential has to be treated in 
perturbation theory (or the wavefunction falls to the origin) Then smearing parameters 
have to be introduced whose meaning is unclear. In contrast, the modern, field-theory 
based formulation cleanly separates the Cornell potential from the transverse exchange. 
This can be suppressed, its range being of order the inverse gluon gap l/E m i n (k), and 
the hyperfine interaction does not need to be so large as to provide ir-p splitting, since 
chiral symmetry accomplishes this |31| and the transverse interaction can be treated in 
perturbation theory. 
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